Scaling properties of the cluster distribution of a critical nonequilibrium model 
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A geometric approach to critical fluctuations of a nonequilibrium model is reported. The two- 
dimensional majority vote model was investigated by Monte Carlo simulations on square lattices 
of various sizes and a detailed scaling analysis of cluster statistical and geometric properties was 
performed. The cluster distribution exponents and fractal dimension were found to be the same 
as those of the (two-dimensional) Ising model. This result, which cannot be derived purely from 
the known bulk critical behaviour, widens our knowledge about the range of validity of the Ising 
00 ■ universality class for nonequilibrium systems. 
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I. INTRODUCTION 



Nonequilibrium models, defined by a stochastic evolution rule, often develop long range correlations in space and 
time. As a consequence of the absence of the microscopic reversibility constraint, critical phenomena in those systems 
exhibit a wider range of universality classes than in the usual thermal phase transitions. A lot of research activity 
' in nonequilibrium systems in recent years was concerned with the characterization of steady states critical behaviour 
CNJ ■ and also with some time-dependent properties (IJ . 

The relation between the singularities in macroscopic quantities and the corresponding internal cluster structure 
in critical equilibrium systems received considerable attention in the seventies and eighties [2-8]. A scaling theory for 
the cluster distribution was presented ||, involving the standard critical indices plus a new exponent describing the 
compactness of the clusters. Geometric clusters were however found to be too large to describe correlations between 
spins. The concept of critical clusters (also called physical clusters or droplets) was introduced || and is the basis of 
fast algorithms of thermal equilibration |J. 

A number of nonequilibrium models have critical behaviour characterized by equilibrium Ising exponents [10-20]; it 
is however not known if their structure on the 'microscopic' level is also Ising-like - as Alonso et al have found for 
a nonequilibrium random field model - or rather a distinct one. The influence of the dynamics on the system steady 
q . state geometric properties is yet to be clarified. 

Cluster-flipping rules have also been used to propose new nonequilibrium models and universality classes [p2[ . 
Nonequilibrium models are defined by a particular, usually local, dynamic rule - modifying this rule often leads 
asymptotically to a distinct steady state. Since local dynamics suffer from critical slowing down, it is desirable to 
devise alternative (faster!) dynamics preserving the statistical properties of a particular steady state. This raises the 
interesting question of how to relate geometric to physical clusters for a particular nonequilibrium model. As a first 
step to answer these questions, we have analysed the cluster structure of a nonequilibrium model that belongs to the 
two-dimensional kinetic Ising universality class, the 2c? majority vote model p^ , pT] , pO]] . The results are compared with 
those obtained in the Ising case. 

II. MODEL AND FORMALISM 

'Spin' variables (voters) {at} — ±1 are defined on the sites of a square lattice and evolve in time by single-spin flips 
with a probability Wi given by 



Wi = i 1 - - 2q)S 



(1) 



where S(x) — sign(x) if x ^ 0, 5(0) = and the sum is over nearest neighbours of <jj; q is the control parameter 
and we will retrict to < q < 0.5 (alignment with neighbours is favoured). At q = q c = 0.075 (for a square lattice), 
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the system undergoes a nonequilibrium phase transition from a ferromagnetic (low q) to a paramagnetic (higher q) 
state; the (static and dynamic) critical behaviour is Ising-like |l(],[ll],^(| . 

In what follows, a cluster is a connected set of spins such that each spin has at least one nearest neighbour in 
the same state. The number of spins in the cluster is the cluster size (or mass) s; spins with at least one neighbour 
with the opposite orientation belong to the cluster perimeter. By analogy with the Ising case, one expects that at 
criticality there are non-spanning clusters of all sizes, i.e. the cluster size distribution n s (average number of clusters 
of size s per lattice site) has power-law form 0] 

n s ~ s- T (2) 

The average perimeter of s-size clusters (per spin) P s scales as 

Ps ~ s a (3) 

where a measures the degree of compactness of the clusters and varies between 1 (ramified clusters) and 1 — \ 
(compact clusters in d dimensions). Corrections to scaling are expected for q slightly away from q c and take the form 

n.(«,e)~a- r /(a»e) (4) 



P s (s,e)~s< 7 g( S ye) (5) 

with e = q "~ q ■ Formulas above are valid for 1 <C s « i" 1 , where L is the system linear size. Assuming expressions 
(2)-(5) hold for the present model, it is straightforward to see that, as in the equilibrium case jH|, the distribution 
exponents t, a and y are related to the usual critical indices by 

/3=— (6) 
V 

a = l- T -^l (7) 

y 

so one additional exponent is sufficient to describe the scaling properties of the non-spanning clusters. For an 
Ising system the spanning cluster is compact in the ferromagnetic phase but becomes a fractal (fractal dimension 
D) at the critical point. In contrast to percolation, where D is related to the euclidean dimension d through critical 
exponents of the same model , for Ising or Potts models it was shown |t],||,^5],|6| that D is independent from the 
corresponding bulk exponents. It can be obtained from the second moment of the (finite) cluster distribution (the 
percolative susceptibility) H 

Xp = Y.^s~e d -™ (8) 

or, for a finite system at criticality, 



Xp (q c ,L)~L- d+2D . (9) 

A direct measurement of the fractal dimension D of the 'infinite' cluster is made by counting the mass within a 
distance r from its center of mass 

s(r) - r D . (10) 
or, alternatively, by using the relation between average size of the 'infinite' cluster Sqo and lattice side L 

Soo~L D . (11) 



Assuming that the finite clusters have the same fractal dimension |j 24 2^,|2^], the dependence of cluster size on 
mean radius of gyration R s also yields D 

s ~ Rf. (12) 
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III. MONTE CARLO SIMULATIONS AND RESULTS 



We have performed simulations of the model described in the previous section on square lattices of 32 2 , 40 2 , 64 2 , 
128 2 and 200 2 sites with periodic boundary conditions for several values of q < q c . After reaching a stationary state 
for each value of q, the system evolution was monitored for up to 1.2* 10 6 MCS; configurations were analysed every 200 
MCS. Clusters of either 'aligned' or 'reversed' spins (compared to the dominant sample orientation) were identified 
and counted separately, together with the respective perimeters. Double logarithmic plots of n s (s) and P s (s) - for 
reversed spins - are displayed in figures 1 and 2, respectively, for several q values; the slopes of the corresponding 
linear regions are consistent with the Ising exponents r = 2.05 and a = 0.68. From fig. 2 we notice that the perimeter 
scaling law (3) is only valid for intermediate size clusters, larger clusters being more ramified than small ones, as 
found in the Ising case |^|. A better estimate of r, together with the correction to scaling exponent y, are obtained 
from the best scaling fit to the data of the function B 



as shown in fig. 3. This data collapse yields r = 2.05 ± 0.005 and y — 0.445 ± 0.01. 

The study of the fractal dimension of finite (inverted) clusters at q c is presented in fig. 4 for the largest system size 
(L = 200). In order to reduce statistical scatter, the results were collected in bins of growing size (bin width increases 
by a factor of 1.5 for each successive interval). The best fit to the data yields D = 1.938 ± 0.003; smaller systems 
produced slightly higher values, so we estimate D = 1.945 ±0.01. This is in agreement with the value D = 1.95 ±0.02 
obtained from the finite size relation (9) for L = 32, 40, 64, 128 and coincides, within error, with the exact Ising value 



Fig. 5 shows a direct evaluation of the fractal dimension of the percolating cluster at q c . The best fit for the L = 200 
system gives D = 1.970 ± 0.003 to be compared with D = 1.968 obtained from relation (11) for L = 64, 128 and 
200. This slightly higher value of D obtained from the infinite clusters data is interpreted as a finite size effect: even 
for the largest lattice, the magnetization is finite at q c and the percolating cluster was allways found to have the 
dominant orientation of the sample. Hence, we were probably not probing the incipient infinite cluster but rather a 
more compact object. 



We have performed a detailed numerical study of the statistical properties of geometric clusters close to the critical 
point of a two-dimensional majority vote model. Scaling was confirmed with distribution exponents r = 2.05 ± 0.005, 
a = 0.68 ± 0.01 and y — 0.445 ± 0.01, indistinguishable from their equilibrium Ising values, as reported in reference 
plf for another model in the same universality class. The fractal dimension associated with the cluster distribution 
at criticality D = 1.945 ± 0.01 was also found to coincide with its Ising value D = A direct determination of D 
from the geometric properties of the 'infinite' cluster produced a slightly higher value, which we believe is due to the 
finite system size. The percolative susceptibility diverges with an exponent ^ p = — 2 + 2D = 1.90 larger than the Ising 
value 7 = 1.75, indicating that, as expected, purely geometric clusters do not coincide with critical ones. Coniglio and 
Klein's rule |^| to reduce geometric to physical clusters - introducing 'active bonds' with a temperature-dependent 
probability - is not easily adaptable to this system. The difficulty to relate the probability for an 'active bond' with 
the noise parameter q is due to the non-hamiltonian nature of the model. The study of the majority vote phase 
transition as a percolative transition of (suitably defined) physical clusters is a challenge for future work. 



We thank J M B Oliveira for his help in the implementation of the cluster counting algorithm and J F Mendes for use- 
full discussions. This work was partially financed by Praxis XXI (Portugal) within project PRAXIS/2/2. l/Fis/299/94; 
M C has benefitted from a junior research grant (BIC) from Praxis XXI. 
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D = 1.947 §. 
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V. FIGURE CAPTIONS 



Fig.l - Double logarithmic plot of the cluster size distribution. From bottom to top q — 0.070,0.071,0.072,0.074 
and 0.0747 (successive curves have been shifted one unit up for clarity). The straight lines have slope —2.05. Lattice 
size L = 200; averages over 5000 or 6000 (q = 0.0747) samples. 

Fig. 2 - Log-log plot of the perimeter length P s as a function of size s. Same conditions as in fig.l. The dotted lines 
have slope 0.68. The inset shows a zoom of the linear region for the five q values; the straight line is a least-square fit 
and has slope 0.683. 

Fig. 3 - Scaling function F versus eS y for r = 2.05 and y — 0.445 (see text). The arrows indicate end of data for 
(left to right) q = 0.0747,0.074,0.072,0.071 and 0.070. Lattice size L = 200. 

Fig. 4 - Cluster size v.s. mean radius of gyration in double logarithmic scale for L = 200 and q = q c ; data binned as 
described in the text. The best fit line has slope 1.938. Averages over 5000 samples. 

Fig. 5 - Mass within distance r from center of mass, s r = s(r), for percolating clusters at q c . Averages over 5000 
samples. System sizes L = 64(D), 128(0) an d 200(+). The straight line has slope D = 1.970. 
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